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Study of 6 dependence in Yang-Mills theories on the lattice 
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We discuss the use of field theoretical techniques in the lattice determination of the free energy dependence on the 8 angle in SU(N) 
Yang-Mills theories. 



1 Introduction 

The dependence of the free energy density F (9) of Yang- 
Mills theories on the 9 angle is the subject of ongoing 
theretical debate. F{9) is defined as: 

exp[-VF(9)] s J [dA] e~^ %m e ieQ (1) 

where = ^F"(x)F^ Y (x) is the usual Yang Mills la- 
grangian and Q = Jd 4 x q(x) is the topological charge, 
with the topological charge density q(x) denned as 

«(*) = ^WrC« F P<rM = <W*) , (2) 

where K^{x) is the Chern current 

K, = (d p A a a - l -gf^A b p A^ (3) 

The determination of F(0) is a typical non-perturbative 
problem and lattice QCD is in principle a natural tool to 
deal with it, but the complex nature of the euclidean action 
for 9^0 forbids the use of standard Monte Carlo simu- 
lations. However many interesting physical aspects can be 
analyzed by studying F(8) for small values of 6. The terms 
in the Taylor expansion of F(0) around 8 = 

CO 

F{8) = F(0) + V -F W (O)0* (4) 

ti kl 

are related to the connected moments of the topological 
charge distribution at 9 = 0: 

J*»(0) S JLf(0| m = -*<^, (5) 

which can be determined by standard lattice Monte Carlo 
simulations 1 . 

1 Other approaches for a determination of F(0) have been tried, based on 
a Fourier transform of the topological charge distribution or on extrapola- 
tions from simulations at imaginary values of 6 [Q] 



The quadratic term is proportional to the topological sus- 
ceptibility, x — (Q 2 )/V> which is expected to be + to the 
leading order in l/N c in order to solve the so-called U(l) 
problem [ |2| [3) . It has been already extensively studied 
by lattice simulations (see Refs. [|4][5] for recent reviews): 
since the study of topological quantities on the lattice is al- 
ways non-trivial, several methods have been used: cooling, 
the field theoretical method, and fermionic methods based 
on the index theorem, all giving consistent results for^ and 
in agreement with the Witten-Veneziano formula relating^ 
to the 7]' mass. 

It has been argued that, for 9 < n, F{9) is almost quadratic 
in 9, with 0(9*) corrections suppressed by powers of 1 /N c [ 
|6|[7). In order to verify this conjecture it is interesting to 
obtain a lattice determination of the quartic term in the ex- 
pansion of F(9) and measure its relative weight with re- 
spect to the quadratic one. This has been done using the 
cooling method [ 8 1 and the field theoretical methd [|9j. In 
the following I will illustrate the details of the lattice deter- 
mination in the framework of the field theoretical method 
and compare the results with those obtained by the cooling 
method. Further details can be found in Ref. [ 9 1. 

2 The method 

On the lattice it is possible to define a discretized gauge 
invariant topological charge density operator qL(x), and a 
related topological charge Ql = 2 1 qi(x) (with the sum 
extended over all lattice points), with the only requirement 
that, in the formal (naive) continuum limit, 

q L (x) a ~°a 4 q(x) + 0(a 6 ) , (6) 
where a is the lattice spacing. A possible definition is 

-1 ±4 

qdx) = X WrTr(lWx)IWx)) ' ( ? ) 

pivpcr-±l 

where Il /JV (x) is the usual plaquette operator in the fiv plane, 
e^ypo- is the standard Levi-Civita tensor for positive direc- 
tions and is otherwise defined by the rule e^pa = -e^- M ) VptT . 
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A proper renormalization must be performed when going 
towards the continuum limit, like for any other regular- 
ized operator. For instance, in spite of the formal limit in 
Eq. 0, the discretized topological charge density renor- 
malizes multiplicatively [ 10 1: 



q L (x) = Z(/3)a 4 (J3)q(x) 



■ 0(a 6 ) 



(8) 



with a multiplicative renormalization constant Z(J3) which 
is a finite functic 
proaching 1 as 



is a finite function of the bare coupling /? = IN/g^, ap- 



Further renormalization constants, relating lattice to con- 
tinuum quantities, may appear when defining correlation 
functions of the topological charge. For example, in the 
case of the topological susceptibility, 



X = 



<e 2 > i ,4 



V 



■I 



(fx (q(x)q(Q)) , 



(9) 



one in general is not guaranteed that the lattice definition 



XL = Yj<^(x)q L (0)) 



(10) 



the background topological content, in order to remove the 
renormalizations, i.e. Z — > 1 and M — > 0). 

Improved (smeared) operators [ 1191 can be used in order 
to reduce the renormalization effects, thus leading to im- 
proved estimates ofx- 

Let us now turn to the problem of the determination of 
the connected quartic moment, (Q 4 ) c = (Q 4 ) - 3(£> 2 ) 2 . 
It is clearly necessary to first understand how the lattice 
expectation value (QT) renormalizes with respect to the 
continuum {Q 4 }. Apart from an obvious Z(J3) multiplica- 
tive renormalization, there will be additive renormaliza- 
tions coming from contact terms appearing in 



<2i> = j ' d 4 xi ■ ■ ■ A 4 x A {q L (xi) . . . q L (x 4 )) 



(13) 



as two or more charge densities come to the same point 
(x,- ~ xj for some i, j). 

It can be shown [ 9 1 that a quite natural assumption for the 
renormalization rule is the following 



(Q 4 L ) = Z(/3) 4 (Q 4 ) + M 4a (J3)(Q 2 ) + M 4 ,o08) 



(14) 



satisfy the correct continuum prescription for the contact 
term arising in Eq. as x — > 0, and this leads to the ap- 
pearance of an additive renormalization, so that the lattice 
quantity xl is related to the continuum quantity^ by 



XL 



Z(J3) 2 a 4 (J3) X + M(J3) 



(ID 



The idea of the field theoretical method for the determi- 
nation of x is to compute and then subtract the renormal- 
ization constants from the lattice quantity Xl> following 
Eq. il It . The determination can be performed numerically 
using the so-called heating method r llOlll llll2lll3U14l[T5l 
U6lll7lfT8l : the renormalization constants are related to the 
UV fluctuations living on the scale of the lattice spacing 
that, close enough to the continuum limit, are effectively 
decoupled from the topological modes living on physical 
scales. It is thus possible to create samples of configura- 
tions of fixed topological background with the UV fluctua- 
tions thermalized, by heating a semiclassical configuration 
of initial charge Q. Expectation values on this samples give 
the necessary information: 



This assumption can be shown to be theoretically sensi- 
ble and allows a straightforward extension of the heating 
method to determine the two new renormalization con- 
stants M 4i2 0S) and M 4 ,o(y6). 

Indeed, the measurement of {Q 4 L ) on a sample of configu- 
rations with background topological charge Q gives 



(Q 4 L ) = Z 4 (P) Q 4 + M 4 , 2 C0) Q 2 + M 4)0 OS) 



(15) 



and if the measurement is repeated in at least two different 
topological sectors, sufficient constraints are obtained to 
determine M^dp) an d M4 0Q3), If more topological sectors 
are used, there is an excess of constraints which can be 
exploited as a non-trivial test of the method. 

In the general case of the n-th order correlation function, a 
natural extension of Eq. (I14> is the following: 



«/2 

(Q'l) = Z"(Q") + M n ^ 2h (Q"- 2h ) 

h=\ 



(16) 



(Ql) 

(Ql) 



Z(J3)Q 

Z(fi) Q 2 + V M(J3) 



(12) 



and its validity can be discussed along the same lines as for 
«=4[|9). 



3 Results 



(this is also, in some sense, the idea at the basis of cooling, 
in which the UV fluctuations are suppressed by a process of 
local minimization of the action, without hopefully altering 



The method has been applied to the case of S U(3) pure 
gauge theory, on a 16 4 lattice at /3 — 6.1, using 1-smeared 
and 2-smeared operators 1 I19II16I . 
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We have collected five different samples of configurations, 
one thermalized in the Q — sector (around the zero 
field configuration), two in the Q — 1 sector (thermal- 
ized around two different semiclassical configurations of 
topological charge one) and two in the Q — 2 sector (ther- 
malized around two different semiclassical configurations 
of topological charge two). The semiclassical configura- 
tions have been obtained by extracting thermalized con- 
figurations with non-trivial topology from the equilibrium 
ensemble at = 6.1 and then minimizing their action 
by a usual cooling technique. All the five samples have 
been obtained by performing about 3000 heating trajecto- 
ries around the semiclassical configurations, each trajec- 
tory consisting of 90 heating steps; 6 straight cooling steps 
have been applied on heated configurations to check that 
their background topological content did not change. 

We have then measured the expectation values (Ql), (Q 4 L ), 
and also (Ql)/Q where Q + 0, over the five samples. We 
have reported the results in Table 1 for the 1 -smeared oper- 
ator and in Table 2 for the 2-smeared operator: expectation 
values obtained on samples with the same Q turned out 
to be equal within errors, as they should, and we have re- 
ported in the tables only their weighted averages. Those 
data have then been used to perform a fit to Eqs. fl!2i 
and Jl 5i obtaining finally the values of the renormalization 
constants reported in Table 3. 



Table 1. Expectation values measured in different topological 
sectors for the 1 -smeared operator. 



Q 


Z = (Q L )/Q 


(Ql) 


(Ql) 







0.311(12) 


0.290(20) 


l 


0.416(6) 


0.4785(60) 


0.630(15) 


2 


0.413(5) 


0.9626(80) 


1.973(50) 



Table 2. Expectation values measured in different topological 
sectors for the 2-smeared operator. 



Q 


Z = {Ql)/Q 


(Ql) 


(Ql) 







0.208(10) 


0.124(10) 


l 


0.544(5) 


0.489(5) 


0.556(12) 


2 


0.542(4) 


1.314(8) 


2.77(6) 



The equilibrium values {Q 2 L ) and {Ql), which are reported 
in Table 4, have been measured on a sample of 300K con- 
figurations separated by five updating cycles, each com- 



3 

posed of a mixture of 4 over-relaxation + 1 heat-bath up- 
dating sweeps; the reported errors have been estimated by 
a standard blocking technique. 

Using Eqs. Jl 11 1 and fl!4i we can compute {Q 2 ) and {Q 4 }, 
obtaining the results reported in Table 4. It is interesting to 
notice that the values obtained for the 1 -smeared operator 
and for the 2-smeared operator are in good agreement, as 
they should, confirming the correctness of the method. 

We can finally determine <g 4 > c = (Q 4 ) - 3(Q 2 ) 2 , obtain- 
ing (Q 4 ) c = 0.32 + 1.80 for the 1-smeared and (Q 4 ) c = 
0.66 + 0.90 for the 2-smeared operator, leading to bi - 
-0.012(62) and b 2 = -0.024(32) for the 1-smeared and 
2-smeared operator respectively, in agreement with the de- 
termination reported in Ref. 



Table 3. Values of the renormalization constants obtained re- 
spectively for the 1-smeared and 2-smeared operators, by using 
the results reported in tables 1 and 2 and performing a best fit to 
Eqs. 03 and {BJ- 



z 


VM 


M 4 ,o 


M 4 ,2 


0.414(4) 
0.543(5) 


0.315(6) 
0.211(5) 


0.336(16) 
0.377(15) 


0.289(16) 
0.124(9) 



Table 4. Expectation values measured at equilibrium and results 
obtained for the renormalized quantities, respectively for the 1- 
smeared and 2-smeared operators. 



(Ql) 


(Ql) 


(Q 2 ) 


(Q 4 ) 


0.7121(38) 
0.8776(60) 


1.548(18) 
2.368(36) 


2.312(72) 
2.262(41) 


16.4+ 1.8 
16.02(72) 



4 More on the renormalization effects 

The renormalization constants Z, M, h „, (m < n) which, for 
a given lattice discretization Qi, appear in Eq. d!6i . are in 
principle independent of each other, or at least no simple 
relation exists among them, unless some further hypoth- 
esis can be done about the nature of the UV fluctuations 
which are responsible for the renormalizations. We will 
propose and test an ansatz which will greatly simplify the 
structure of the renormalization constants and will lead to 
a renormalization formula which directly involves the con- 
nected correlation functions, thus allowing a more precise 
determination of b^. 



4 
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An hypothesis about the nature of the UV fluctuations has 
been done in Refs. HI 1 1 1 1 31 . where it was assumed that the 
discretized topological charge density can be expressed as 



q L (x) st [Z + ((x)]q(x) + T](x) 



(17) 



where q(x) is a background topological charge density 
which is determined by physical fluctuations on the scale 
of the correlation length whereas £(x) and j](x) are ran- 
dom variables with zero averages which are determined by 
the short range UV fluctuations and, at least in the con- 
tinuum limit, are expected to be decoupled from q(x), i.e. 
(((x)q(x)) = (r](x)q(x)) = 0. Summing Eq. dl7> over all 
lattice points, the following relation follows for the lattice 
topological charge Qjj. 



Q L =ZQ + Y j ax)q(x) + J 1 . 



(18) 



where rj = Ytx rj(x). We now make the further assump- 
tion that the term £ A Z( x )<l( x ) m Eq. J 1 8t > can be neglected, 
configuration by configuration. This is not unreasonable, 
in view of the fact that q(x) and g(x) are decoupled from 
each other. We will thus assume that 



Ql = Z Q + r] . 



(19) 



where 77 is a random noise with zero average which is 
stochastically independent of Q. 

This assumption has relevant consequences for the struc- 
ture of the renormalization constants. Indeed, using the hy- 
pothesis that Q and r\ are stochastically independent vari- 
ables and that they are both evenly distributed around zero, 
it is easy to verify that the general renormalization formula 
holds: 



n/2 



h=0 



2h 



\Z"- 2h (ff'- zn }(if n ) , 



2h\/lh\ 



(20) 



so that the renormalization relation for (Q" L ) is described 
only in terms of Z and of the correlation functions of the 
noise rj. In particular we have M n ,m — (l)Z m (if~ m ), a rela- 
tion that should be verified on numerical data if our ansatz 
in Eq. dl9> is correct. From the data in Table 3 it can be 
checked that indeed M 4>2 = 6Z 2 (t] 2 ) = 6Z 2 VM, but we will 
now proceed further and check the validity of Eq. J20i up to 
n = 6. The correlation functions of 77 can be determined by 
the heating method using the analogous of Eq. i!5i . which 
up to n — 6 reads: 



(Q\) = Z 2 Q 2 + {i 1 2 ) 

(Qi) = z 4 Q 4 + 6z 2 eV> + (n 4 ) 

(Ql) = Z 6 Q 6 + 15Z 4 Q 4 (7 ] 2 ) + 
+ \5Z 2 Q 2 {n 4 ) + {if) . 



Using the values for (Q 2 L ), (Qi) and (Q 6 L ) obtained in the 
sectors with Q = 0, 1,2 we have performed a best fit to 
Eqs. fl2 li . obtaining the best fit values reported in Table 
5 with ^ 2 /d.o.f. =2 0.34 for the 1 -smeared operator and 
X 2 /d.o.f. =2 0.23 for 2-smeared operator. The fact that the 
values for (Q 2 L ), (Q 4 L ) and (Q b L ) obtained in the various sec- 
tors can be fitted by the simple relations in Eq. 11211 is a 
confirmation of the validity of the ansatz in Eq. d!9i . 

Assuming that Eq. d!9l > is valid, it is possible to write a 
renormalization relation which involves directly the con- 
nected correlation functions. Indeed, it is a general rule that 
the connected correlation functions of a stochastic variable 
(Ql in our case), which is the sum of two variables which 
are stochastically independent of each other (ZQ and r\ in 
our case), are the sum of the corresponding connected cor- 
relation functions, i.e. 



(Q" L )c = Z"(Q") C + (rf) c 



(22) 



Therefore in order to compute (Q") c we need to know, 
apart from Z, only one renormalization constant, (77" ) c , 
which can be easily measured by computing {Q" L ) C on the 
sample of configurations in the Q — sector. (Q") c is then 
given by 



(Q")c 



(Ql)c - (rf)c 
Z" 



(23) 



where (Q" L ) C is measured on the ensemble of configurations 
at equilibrium. 

We have computed (Q 4 L ) C - (Q A j) ~ 3(Q?) 2 on our equilib- 
rium configurations atyS = 6.1, obtaining (Q 4 L ) C = 0.026(7) 
for the 1-smeared operator and (Q 4 L ) C = 0.057(13) for the 
2-smeared operator. We have then computed (Q 4 L ) C on our 
sample of configurations thermalized in the Q = sector at 
P = 6.1 , obtaining (i] 4 ) c = (rj 4 ) - 3(tj 2 ) 2 = 0.006(4) for the 
1-smeared operator and (?7 4 )c = 0.001(2) for the 2-smeared 
operator. In both cases (equilibrium and 2 = 0) errors have 
been estimated by standard jackknife techniques. 

By using Eq. (I23> and the values for Z and (Q 2 ) pre- 
viously obtained, we have obtained (Q 4 ) c = 0.68(24), 
/>2 = -0.024(9) for the 1-smeared operator and (Q 4 ) c 
O.661 15 ), b 2 = -0.024(6) for the 2-smeared operator. 

By making use of the ansatz in Eq. dl9l we have thus made 
determinations which are much more precise than those 
obtained in Section 3. The reason is that Eq. i22\ allows to 
relate (Q") c directly to the connected correlation functions 
of the discretized lattice topological charge, with only two 
renormalization constants involved: this greatly simplifies 
computations and error propagation, thus leading to im- 
proved estimates. We notice that most of the error in the 
final determination of (Q 4 ) c and bi comes from the deter- 
mination of (Q 4 L ) C at equilibrium, which is also the most 
(21) expensive part of the computation in terms of CPU time. 
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The renormalization procedure is thus completely under 
control and numerically non expensive. 

We have also made a determination of bo at /3 = 6.0, again 
on a 16 4 lattice. On a sample of about 300K configura- 
tions obtained at equilibrium and using the same algorithm 
as for /3 = 6.1 we have obtained, for the 2-smeared op- 
erator, (Q 2 L ) = 1.377(7), (Q 4 L ) C = 0.052(23). On a sam- 
ple of configurations thermalized in the Q — topological 
sector by performing about 3000 heating trajectories, each 
composed of 90 heating steps, we have obtained, for the 2- 
smeared operator, (tj 2 ) = 0.308(10) and <77 4 ) c = 0.002(3). 
From these data, using the value Z(J3 = 6.0) = 0.51(2) 
reported for the 2-smeared operator in Ref. r i!6l . we ob- 
tain Z?2 = -0.015(8), which is consistent with the value 
obtained at/? = 6.1. 

Let us close noticing that the value obtained for (?7 4 ) c is 
very small and compatible with zero for both the 1 -smeared 
and the 2-smeared operator. We have also measured (rj 6 ) c 
on the sample of configurations at Q — obtaining (t] 6 ) c = 
0.001(8) for the 1-smeared and {if) c = 0.0005(14) for 2- 
smeared operator (J3 = 6.1), so that 77 behaves with a good 
approximation as a pure gaussian noise. 



Table 5. Values of the renormalization constants, respectively for 
the 1-smeared and 2-smeared operators, obtained by performing 
a best fit to Eq. ETJ. 



z 


(t] 2 ) 


<rj 4 ) 


in 6 ) 


0.415(4) 
0.542(4) 


0.315(6) 
0.211(5) 


0.298(11) 
0.129(6) 


0.462(37) 
0.131(17) 



5 Conclusions 

We have discussed the extension of the field theoretical 
method, already used for the lattice determination of the 
topological susceptibility, to the computation of further 
terms in the expansion of the ground state energy F{8) 
around 8 = 0. 

We have presented numerical results regarding SU(3) pure 
gauge theory, providing a determination of the fourth order 
term in the expansion of F(ff) around 8 = 0. Our determi- 
nation is in agreement with that obtained in Ref. [ 8 1 via the 
cooling method, and confirms that fourth order corrections 
to the simple 8 2 behaviour of F(8) around 8 = are small 
already for N c = 3. 
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